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ABSTRACT 



Context. The determination of pulsation velocities from observed spectra of Cepheids is needed for the Baade-Wesselink calibration 
of these primary distance markers. 

Aims. The applicability of the Fourier-disentangling technique for the determination of pulsation velocities of Cepheids and other 
pulsating stars is studied. 

Methods. The KOREL-code was modified to enable fitting of free parameters of a prescribed line-profile broadening function corre- 
sponding to the radial pulsations of the stellar atmosphere. It was applied to spectra of S Cep in the H-alpha region observed with the 
Ondrejov 2-m telescope. 

Results. The telluric lines were removed using template-constrained disentangling, phase-locked variations of line-strengths were 
measured and the curves of pulsational velocities obtained for several spectral lines. It is shown that the amplitude and phase of the 
velocities and line-strength variations depend on the depth of line formation and the excitation potential. 
Conclusions. The disentangling of pulsations in the Cepheid spectra may be used for distance determination. 

Key words. Line: profiles ~ Techniques: spectroscopic ~ Stars: variables: Cepheids 



1. Introduction 

The method of Fourier disentangling of spectra was developed 
by Hadrava (1995) from the method of cross-correlation (cf., 
e.g.. Hill 1993) to decompose, from a series of observed spectra 
of multiple stars, contributions of the individual components and 
to simultaneously find orbital parameters, or, more generally, to 
fit free parameters of the physics governing the Doppler shifts 
and line-profile variations. One of the advantages of disentan- 
gling compared to cross-correlation is that it does not require a 
template spectrum of a star with similar spectral type or a model 
atmosphere. 

Zucker and Mazeh (2006) introduced their method 
of Template Independent RAdial-VELocity measurement 
(TIRAVEL) for single-lined spectroscopic binaries (SBl). Their 
method uses each exposure from a series of observations as a 
template for cross-correlation with all other exposures. Zucker 
and Mazeh mentioned that the KOREL-code (cf., e.g., Hadrava 
2004b) for Fourier disentangling is also template-independent. 
Nevertheless, they advocated use of their TIRAVEL based on 
the alleged advantage that it assumes the individual radial ve- 
locities to be free variables not bounded by an orbital motion. 
However, KOREL enables the convergence of individual radial 
velocities of any component also. Moreover, KOREL can do this 
not only for SBl, which is an extremely simple case, but also for 
two or more components. The option of free velocities is rarely 
used, because in practice it is more advantageous to take the or- 
bital motion of multiple stars into account and to solve directly 
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for the orbital parameters. Zucker and Mazeh mentioned the case 
of a third component as an example of when free radial veloci- 
ties are needed. However, just in this case the solving for orbital 
parameters of both the close and the wide orbit (taking into ac- 
count also the light-time effect) simultaneously with the determi- 
nation of radial velocities is advantageous, because it checks the 
consistency of the possibly small perturbation with the source 
spectra better than a two-step procedure of determination and 
subsequent solution of the radial-velocity curve. 

One reason to solve in some cases for individual radial veloc- 
ities independent of orbital motion is technical, e.g. in the case 
of an unreliable wavelength scale of the observed spectra, as it 
has been done in the study by Yan et al. (2008). However, a more 
important reason is seen in cases when the observed wavelength 
shifts of spectral lines are not due to overall orbital motion of the 
star but due to some other effects. One such case is the Doppler 
shift caused by pulsations of the stars. The use of Fourier disen- 
tangling in the spectroscopic studies of pulsating stars has been 
outlined by Hadrava (2004a,b). Here we shall demonstrate this 
method in practice in the case of the star 6 Cep. 

The study of either radially or non-radially pulsating stars is 
important as a clue to probe the inner structure of stars. In addi- 
tion, the period-luminosity (PL) relation of Cepheids and some 
other radial pulsators is used as one of the few primary meth- 
ods for distance determination. One needs to calibrate this PL- 
relation, and that can be done by the Baade-Wesselink method 
(Baade 1926, Wesselink 1946). In this method, the spectroscop- 
ically measured pulsation velocity is integrated over the course 
of the period to yield the changes of stellar radius, which, com- 
bined with photometric or interferometric variations, may reveal 



2 



P. Hadrava et al.: Notes on disentangling of spectra 



the distance to the star. However, the spectroscopic measurement 
is compHcated by the hne-profile variations caused not only by 
the radial motion of the stellar atmosphere (which must be inte- 
grated over the visible part of the stellar disc projected locally 
to the line of sight) but also due to changes of other physical 
conditions (e.g. the temperature and density) in line-forming re- 
gions of different lines. Consequently, the radial velocities 
measured by different methods reflect only indirectly the instan- 
taneous pulsation velocity Vp of the stellar surface. It is thus used 
to introduce the so called projection factor p = Vp/vr, which can 
be estimated either theoretically (e.g. Nardetto et al. 2004) or ob- 
servationally (Nardetto et al. 2008 and citations therein). An al- 
ternative method directly matching the Doppler shifts and asym- 
metries of spectral lines with a proper model of line-profiles has 
been introduced by Gray and Stevenson (2007). Their approach 
is, in principle, equivalent to the above mentioned disentangling 
of pulsations (Hadrava 2004a,b). 

We summarize the method of pulsation disentangling and its 
use in Sect. |2] and we test it on our observations of the 6 Cep in 
Sect. [3] In Sect.|4]we discuss possibilities of further development 
of the method. 



2. Disentangling of radial pulsations 

The line-profile variations (LPVs) caused by motion of the 
stellar atmosphere are often modelled by integrating Doppler 
shifted spectra of unperturbed (mostly plane-parallel) model at- 
mospheres over the stellar surface. This approximation is used 
for study of the rotational broadening as well as LPVs due to 
radial or non-radial pulsations even though it is obvious that the 
motion may change the structure of the atmosphere and hence 
also the radiative transfer and the line formation in them. In the 
case of Cepheids, it is observed that different lines formed in 
different layers of the atmosphere have slightly different phase- 
dependence of LPVs (cf. Breitfellner and Gillet 1993, Butler 
1993 etc.), which reflects (besides other effects) the changes of 
velocity gradients within the atmosphere. The basic idea of the 
Baade-Wesselink method also assumes that the spectroscopic 
and photometric or interferometric variations are caused by pul- 
sations of the same stellar surface. However, the dilution of the 
atmosphere and its temperature variations and the presence of 
stellar winds in these stars implies that the motion of layers with 
given optical depths may differ from the local velocity of the gas, 
which then influences the profiles and shifts of the spectral lines. 

A safe way (outlined and followed, e.g., by A. Fokin, 1991, 
2003 etc.) to treat all these effects properly would be to con- 
struct physically self-consistent models of atmospheres of pul- 
sating stars, to calculate the observable quantities (spectra, light- 
curves, visibility functions) and to match the values of free pa- 
rameters (including the distance) to the real data. However, the 
modelling itself is still computationally very demanding, not to 
speak of the inverse problem of solving for the free parameters. 
It is thus worth simultaneously following an alternative way of 
fitting the observations by simplified models, which take into ac- 
count the most important effects and could be modified if a dis- 
crepancy with respect to self-consistent models or with respect 
to observations were found. 

Let us assume now that the observed spectrum I(x, t) (as a 
function of the logarithmic wavelength x and time t) is given by 
the integral 



X 



over the visible part (i.e. where the directional cosine ju > 0) of 
the stellar surface s. If the specific intensity (in the rest frame 
of the moving atmosphere) is independent of s and we expand it 
into a power series of ^ 



I{x, s,iJ, t) I^x, t)^/ , 

k 

the spectrum is given by 
I{x, = ^ f) * aHx, f) , 

k 

where 



A*(x,f) = J/ 



/+'5(jc-y(s,f))'i^s . 



(2) 



(3) 



(4) 



For purely radial pulsations synchronous on the whole stellar 
surface, the local radial velocity v{s, t) - i-iVp{t) is the projection 
of the instantaneous speed Vp of the pulsation. The broadening 
functions for individual modes of the limb darkening then read 



A*(x, f) 



2nR^ 

,1+2 



(5) 



where R - R(t) is the instantaneous radius of the star. The 
Doppler shift which can be measured by the first moment is 
given by the mean value of x of the broadening and it reads 



/ xA''{x, t)dx k + 2 
J A*c/x 



k + 3 



(6) 



In agreement with Getting (1934), the projection factor p is 
thus I for lines without limb darkening and | with linear limb 
darkening equal to 1 . For higher order terms of the series given 
by Eq. (|2]l, the projection factor decreases toward the limiting 
value 1 , which corresponds to the broadening function 



A'^{x,t)~6(x-vp(t)) , 



(7) 



I(x, t) - I yu/(x, s,n, f) * 6{x - v{s, t))d s 



(1) 



i.e. to an absorption occurring in the centre of the stellar disc 
only. This calculation is valid exactly for radial-velocity mea- 
surement by the moment method only. For other methods, like 
the deepest point of the profile, bisectors, the fit by two semi- 
Gaussian curves or the cross-correlation method, the result of 
radial-velocity measurement of the asymmetric lines is less cer- 
tain. It also depends on the instrumental broadening. Because 
the standard Fourier disentangling (which assumes the broad- 
ening function given by Eq. Q only) is a handy method of 
radial-velocity measurement, it is worth investigating its prop- 
erties (and its p-factor in particular) when applied to Cepheids, 
or modifying it by taking into account the proper broadening 
function so that it will avoid the /^-factor completely and will 
directly provide the pulsational velocities. 

It is generally understood that limb darkening is crucial for 
the asymmetry of the lines and thus also for setting the projec- 
tion factor However, determination of its proper value is still 
a problem (cf. Montanes Rodriguez and Jeffery 2001, Marengo 
et al. 2002). The limb-darkening corresponding to the radiation 
in continuum is sometimes accepted also in the lines, arguing 
that it "varies slowly with the wavelength and can be taken to be 
constant over the span of a spectral line" (Gray 2005, p. 436). 
However, the spectral flux in continuum also varies slowly with 
the wavelength, but it cannot be taken to be constant over any 
line because it is just its fast variation across the line-profile that 
is seen as the spectral lines. A simple model of radiative transfer 
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in stellar atmospheres reveals that weak lines are usually domi- 
nated by the central part of the visible stellar disc and hence the 
differences between the radiation in the lines and in continuum 
has the limb darkening close to 1.0 (cf. Hadrava 1997). This is 
confirmed by the observed behaviour of the line strengths during 
eclipses and it is also consistent with the values p ^ 1.3 of the 
projection factor in Cepheids. However, detailed non-LTE mod- 
els of stellar atmospheres show that the limb darkening within 
strong lines (e.g. Balmer lines) is more complex (cf. Hadrava 
and Kubat 2003). It can be represented as a superposition of 
a part without limb darkening, a part with linear darkening, as 
well as of higher order terms, which also may include an emis- 
sion (especially at the outer edges of the stellar disc if the spheri- 
cal instead of the plane-parallel symmetry is taken into account). 
We shall thus investigate separately the cases of simulated pul- 
sationally broadened lines without any limb darkening {k - 0) 
and with linearly proportional broadening {k - 1). We decon- 
volve these profiles (and in the next Section also real observed 
spectra) by disentangling with broadening functions for the same 
cases ik - 0, 1) and by the standard disentangling (i.e. k - oo). 

The generalization of disentangling for line-profile variabil- 
ity has been described by Hadrava (1997 and 1998 for the line- 
strength variability and the changes in shape of line-profile, 
respectively; unfortunately, the formulae were incorrectly pro- 
cessed in the print of the later paper). The explicit form of the 
pulsational broadening and its Fourier transform was given for 
the case of limb-darkening equal to 1 by Hadrava (2004a,b). The 
Fourier transform (from the space of functions of the variable x 
to functions of y) of this profile (Eq. (|5]) for ^ = 1) has the form 

^' ^) = ^ITT fexp(/3'Vp)(2 - liyvp - y^v]) -l] . (8) 

Here we need also the case with zero limb-darkening (k - 0), 
for which the Fourier transform reads 

AV, t) = 4^ fexp(i>p)(l - iyvp) - ll , (9) 

and also the standard disentangling, for which K'^iy, f) ^ 
exp(!yvp). The changes of R(t) in these formulae influence the 
flux in the lines in absolute units proportionally also to the vari- 
ations of the flux in continuum. Hence, they should not affect 
the line-strengths in the rectified spectrum (i.e. normalized to 
the continuum), unless there is also light from a binary compan- 
ion or a background star present in the spectrum. However, be- 
cause the changes of the effective temperature and other param- 
eters of the stellar atmosphere in the course of the pulsational 
cycle generally affect the line-strengths, we substitute the term 
2nR^ = (k + 2)s(f), where s(t) is a multiplicative line-strength 
factor given by the integral of Eq. ^ over x (which scales the 
broadening normalized to a unit integral). 

We have calculated several sets of simulated data by di- 
rect integration in the jc-space to always convolve one fixed 
Lorentzian line-profile {<p(x) ~ ((x - xq)^ + y^)"') with broaden- 
ing functions either A" or A' for instantaneous values of velocity 
Vp{t) of a harmonic (i.e. sinusoidal) pulsation at 20 values of t 
uniformly covering the period of the pulsation. The Lorentzian 
profile was chosen because it contains both a narrow core and 
wide wings. An example of such Lorentzian profiles with the in- 
trinsic semi-halfwidth y corresponding to 30km/s broadened by 
a" with a semi-amplitude of the pulsational velocity A',,,,, equal 
to lOOkm/s is presented in Fig.[T] (The step of lOOkm/s is marked 
by the ticks on the jc-axis.) It can be seen that the profiles (cf. the 



/ 















































-100 \ 





Fig. 1. Standard disentangling (with A^^) of simulated 
Lorentzian profiles broadened by radial pulsations (with A°.^, 
q = 0.3). The simulated profiles are drawn by thick lines (blue in 
the electronic version), their fits using disentangling by thin (red) 
lines and the mean disentangled profile by the bottom (green) 
line. 



thick lines with offsets proportional to the phase) have the high- 
est asymmetry at extremes of the pulsational velocity, the high- 
velocity wing always being steeper than the wider low-velocity 
wings reaching the rest wavelength of the line. Qualitatively the 
same feature also results if the pulsational broadening is applied 
to non-Lorentzian line profiles and is in agreement with the line- 
profile variations observed in Cepheids (cf., e.g., Breitfellner and 
Giflet 1993, Nardetto et al. 2006, Gray and Stevenson 2007). It is 
used to characterize such profiles quantitatively by bi-Gaussian 
fits (cf. Nardetto et al. 2006), which, however, are not physically 
substantiated. 

Figure[T]also shows the disentangling of these simulated pro- 
files by the standard KOREL disentangling (i.e. A°°). For spectra 
containing lines of two or more component stars, the disentan- 
gling simultaneously decomposes them and fits the orbital pa- 
rameters, while in the present case it leads to a simpler problem 
of fitting all the input spectra by a mean profile (which is drawn 
by the lowermost line) scaled in the line-strength and shifted in 
X. The best fits to the individual input profiles are overplotted in 
Fig. [Ilby the thin lines. The disentangled mean profile is sym- 
metric here (which need not be exactly the case if the pulsational 
velocity includes some higher overtones in addition to the basic 
sinusoidal mode, or if the spectra do not cover the period uni- 
formly). Yet, it reproduces the input profiles relatively well, so 
that for real data the difference between the true profiles and the 
fit by a simplified model may be hidden in the noise. The co- 
incidence is even better for input profiles broadened by A' and 
disentangled by A*". In both cases the coincidence is worse if 
the amplitude of the pulsational velocity exceeds more the in- 
trinsic width of the line. The disentangled profile is wider than 
the intrinsic profile (this is well seen e.g. in the uppermost in- 
put line, for which Vp - 0), whenever the disentangling is per- 
formed using a smaller broadening (i.e. A* with higher k) than 
the broadening used for the simulation of the data. If the input 
profiles are disentangled with the same broadening for which 
they were simulated, their fit as well as the reconstruction of the 
intrinsic profile is perfect within the numerical precision of the 
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Table 1. Projection factors for simulated line-profiles 



A" A'. 

Sim Sim 

q A*;. A',. A" A" A' A7 

ais ctis ais dis ais ats 

0.1 0.9999 1.0718 1.3942 0.9426 0.9999 1.2556 

0.3 0.9999 1.1074 1.4624 0.9100 0.9999 1.3064 

1.0 1.0001 1.1229 1.4945 0.8918 1.0001 1.3297 

3.0 1.0013 1.1262 1.5012 0.8905 1.0015 1.3348 

10. 1.0006 1.1253 1.5007 0.8895 1.0003 1.3341 



simulation. The disentangled line-strength factors are higher (up 
to nearly -1-0.08) for phases with low pulsational velocities and 
smaller (nearly -0.08) at extremes of the velocity. 

The disentangling of this and other simulated datasets is 
performed here with velocities bound to a sinusoidal pulsation, 
which was also assumed in the creation of the data. This is done 
by formally assuming that the velocity obeys a circular motion. 
The ratio of the true amplitude Ksi„, of the velocity chosen for the 
simulation and its disentangled value K^hs gives a mean value of 
the projection factor p for a particular combination of the broad- 
ening functions. The values of pulsational velocities which give 
the best fit of the mean disentangled profile with the individual 
input profiles are also calculated and they reveal that the projec- 
tion factor is slightly phase-dependent. However, the differences 
are very small - in the particular case shown in Fig. [T] the free 
radial velocity is about 0.62km/s (i.e. nearly 1%) higher than 
the harmonic one at their extremes and about 0.47km/s lower at 
medium (non-zero) values of velocity. The profiles reconstructed 
with both velocities are drawn in Fig.[T]by the thin lines, but they 
are not distinguishable within the precision of the graphics. We 
thus neglect these differences which are of the order of the pre- 
cision of the computation and we give the mean values of pro- 
jection factors p - KsimlKdis only for several combinations of 
broadening functions used for the simulation (A^^^^j, ^ = 0, 1) and 
the disentangling (A*.^, k - 0, 1, oo) and the velocity- amplitude 
to line- width ratios q - y/Ksm in Table[T] 

It can be seen from these results that for lines broader than 
their Doppler shift {q - ylKsim > 1) the standard disentangling 
(with A^^) yields /^-factors agreeing within the numerical errors 
(which are of the order 10"^ here) with the moment method of 
the radial-velocity measurement, i.e. | or | for lines with unit 
(A'. ) or zero (A°. ) limb darkening. For lines with intrinsic 

^ Sim' ^ Sim' ^ 

widths y smaller than the pulsationally induced shifts and asym- 
metries, the standard disentangling is more sensitive to the posi- 
tion of the deeper parts of the profile and the /^-factor decreases 
slightly closer to the value 1 . The disentangling with the proper 
limb-darkening in the line (i.e. AV forA'. andA°. for A". )has 

^ ^ ais sun ais sim' 

/:> = 1, it means that it directly provides the pulsational velocities 
and it is desirable to use it for the Baade-Wesselink calibration. 
The use of improper broadening functions (i.e. A'^.^ for A'j.^^^ or 
vice versa), however, results in an error of about 10% in radial 
velocity. We thus need either to find the proper limb darkening 
across the line-profiles from detailed model atmospheres, or to 
distinguish which model fits the observed line-profiles better. 

3. Disentangling of the observed spectra 

To test the disentangling of pulsations on real data we started 
spectroscopic observations of 5 Cep using the 700-mm camera 
of the spectrograph in the Coude focus of the Ondfejov 2-m tele- 
scope equipped with LN2-cooled SITe CCD detector ST-005A 



(2030 X 800 15-jum pixels). Sixty nine medium-resolution spec- 
tra (7? ~ 13 000) with a linear dispersion of 17 A/mm (0.25A/pix) 
in Ha region (6250-6770 A) obtained between August 19, 2008 
and April 16, 2009 (mostly by M. Slechta) are used in this study. 
See on-line Table[3]or|4]for the journal of observations. 

The spectra were reduced in using the standard pack- 

ages ccdproc, doslit and rv (for more details of the process- 
ing see Skoda and Slechta, 2002). 

The spectral region around the Ha-line contains many at- 
mospheric water-vapour lines which complicate the measure- 
ment of stellar spectra by standard methods (cf., e.g.. Kiss and 
Vinko 2000). However, Fourier disentangling with variable line 
strengths is not only suitable to remove the telluric lines (cf. 
Hadrava 1997, 2004a,b, 2006a), but at the same time it en- 
ables one to use them for an additional check or correction of 
the wavelength scale similarly to their use in classical methods 
(cf., e.g., Butler and Bell 1997). As the first step of the disen- 
tangling we chose the spectral region 6511-6521 A sampled in 
1024 bins (i.e. with a step in radial velocity of 0.45km/s per bin) 
to find the line-strength coefficients for the telluric lines. The 
Doppler shifts of the telluric lines with respect to the heliocen- 
tric wavelength-scale were calculated in the Keplerian approxi- 
mation of the annual motion which is provided by the PREKOR- 
code from the coordinates of the target star (cf. Hadrava 2004b) 
and the telluric lines were disentangled by the standard disen- 
tangling using the broadening function A^^. The stellar lines in 
this region were disentangled as a superposition of two systems 
of lines using A^'^.^ with free radial velocities and each one with 
its own free line-strength factors. To avoid uncertainties in low 
Fourier modes, which could cause anticorrelated distortions of 
the stellar and telluric component continua, we used disentan- 
gling constrained by a template (cf. Hadrava 2006b) for the tel- 
luric lines. The template had been calculated by disentangling 
spectra of the star 68u Her (i.e. spectroscopic binary HD 156633, 
also taken with the Ondfejov 2m-telescope), which gives a tel- 
luric spectrum with a satisfactorily flat continuum. The differ- 
ences between the prescribed annual motion of the telluric lines 
and the radial velocities of these lines disentangled as the best 
fits of individual exposures were then used as corrections of the 
wavelength scale in the preparation of spectral regions for subse- 
quent disentangling of the stellar lines. Owing to application of a 
method of enhanced precision (Hadrava 2009), these corrections 
were found with sub-pixel resolution. These wavelength correc- 
tions can be applied to exposures with sufficiently strong telluric 
lines only. In our case, the depth of telluric lines in exposures 
where they are weakest is about one half of their mean depth, so 
that the correction could be applied to all exposures. 

To disentangle the pulsational velocities in the observed 
spectra of 6 Cep we chose first several narrow spectral regions, 
each one containing a single dominant spectral line listed in 
Table |2l These regions were sampled in 256 bins each with a 
step of radial velocity per bin given in the Table {rvpb in km/s). 
In some of the regions, blends with some weak lines can be seen, 
which partly decrease the quality of the disentangling, however, 
their influence can be neglected. The relatively low spectral reso- 
lution of our original data (about 12 km/s) provides only a rough 
sampling of line profiles in individual exposures, which have 
half-depth widths comparable to the amplitude of the Doppler 
shifts, i.e. only about 4 to 6 times larger. The disentangled line 

' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA), Inc., under cooperative agreement 
with the National Science Foundation 
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A [A] rvpb multiplet ^ + e [eV] EW [A] 



6265.14 
6343.71 
6347.10 
6355.04 
6358.69 
6400.01 
6456.39 
6663.45 



0.7 
0.8 
0.7 
0.6 
0.6 
1.5 
0.8 
1.0 



Fe I (62) 
Cal(53) 
Si II (2) 
Fe I (342) 
Fel(13) 
Fel(816) 
Fe II (74) 
Fel(lll) 



2.18 
4.44 

8.15 + 8.12 
2.83 
0.86 
3.60 

7.87 + 3.90 
2.42 



0.107 
0.184 
0.245 
0.110 
0.134 
0.263 
0.443 
0.188 



19.89 
20.40 
18.88 
19.82 
18.28 
21.42 
17.78 
20.67 



0.402 
0.531 

-0.203 
0.508 
0.712 
0.319 

-0.069 
0.384 




Fig. 2. Phase-dependence of pulsational velocity v in km/s for 
individual spectral lines disentangled with A^.^. The phase (f> is 
labeled in cycles (i.e. the non-integer part of the ephemeris E). 
Lines of neutral atoms are marked by open circles (blue in elec- 
tronic version), lines of ionized atoms by (red) crosses. 



profiles are significantly smoother due to the averaging of a large 
number of the exposures. Nevertheless, the limited quality of 
the data does not enable us to convincingly decide which of the 



models A'l. 

dis 



Al,. and best fits the observed line-profile vari- 



ations or even to search for their best linear combination, which 
could be expected for more general limb-darkening within the 
line-profile. We thus performed the disentangling of all regions 
for each of these models separately. In all regions, the integrated 
(O-C)^ of the fitted spectra was the largest for the model A^^, 
while the residual noise for the pulsational broadenings A°.^ and 
A'^.^ was nearly the same (within about 4%) without any evident 
regular preference of one model or the other. We thus illustrate 
the results on the case A^.^, which corresponds to the unit limb 
darkening and is thus the most advantageous from the theoretical 
point of view (see above). 

Figure |2] shows the pulsational velocity disentangled using 
the model A\. which is given in the on-line Table [3] For the 



Fig. 3. Phase-dependence of line-strength factors s for individ- 
ual spectral lines disentangled with A^.^. (The meaning of sym- 
bols is the same as in Fig.|2l) 



purpose of visualizing the results (i.e. calculation of phase) we 
use the ephemeris 



HJD = 2454697.0009 + 5.366341 x E 



(10) 



Our data do not allow us to check or improve the period, and the 
scatter in its published values (e.g. P = 5'?3663159 by Moffett 
and Barnes, 1985, or 5'?3662351 according to the quadratic 
ephemeris by Berdnikov and Ignatova, 2000) is not substantial 
for our present purpose. The reference epoch is chosen to closely 
precede our observations and coincide with the second exposure, 
which is relatively close to the minimum radial velocity. 

The fit of the observed spectra by the standard disentangling 
(A^^) is invariant with respect to adding a constant to all dis- 
entangled radial velocities and a simultaneous shift of the dis- 
entangled spectrum in the logarithmic wavelength for the corre- 
sponding value. If the position of the stellar spectrum disentan- 
gled from all exposures (transformed first into heliocentric wave- 
length scale) is fixed to the laboratory wavelengths of identified 
stellar lines, the measured radial velocities correspond to instan- 
taneous heliocentric radial velocities. Neglecting the systematic 
errors caused by the line-profile distortions, the heliocentric ra- 
dial velocity of the star could be estimated as the velocity aver- 
aged over the pulsational period 



PJo 



v(f)dt 



(11) 



Having the values of velocities measured with some er- 
rors in several discrete times only, we fit them using least- 
squares method in the standard manner by Fourier series (cf. 
Schaltenbrand and Tammann 1971, or Moffett and Barnes 1985, 
etc.) 



v(f) 



v'o + ^ Vk sindnkt/P + (pt) ■ 



(12) 



k=l 



The mean velocity v is then given by the first term vq. We have 
chosen K - 5, because the residual (O - C)^ of the fit decreased 
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significantly with each additional harmonic term up to this value, 
but its value remains practically constant (given by the rms. error 
of the measurements) for higher degree of the expansion. We 
also used the same Fourier series for the line-strength factors. 

Unlike the standard disentangling, the disentangling with 
broadening functions A'J. and Al. is sensitive to the value of 

dis dis 

the true radial velocity of the star's centre of mass. At the phases 
when the pulsational velocity is equal to zero, the consequent 
line-profile distortion disappears, the instantaneous line-profile 
coincides (within the errors) with the disentangled one and its 
Doppler shift is given by the overall motion of the star only. The 
Doppler shift of the disentangled spectrum should thus corre- 
spond to the intrinsic radial velocity of the star However, be- 
cause a systematic shift of all pulsational velocities can be rela- 
tively well compensated by a shift of the disentangled spectrum, 
the convergence of the solution to the true radial velocity is very 
slow and the value of the velocity is poorly defined, in partic- 
ular for the rough sampling of the line profiles in our spectra. 
(To speed up the convergence we introduced an additive term in 
the velocity, which is possible to converge explicitly by the sim- 
plex method.) The values for which we achieve the best fit of 
the observed spectra differ somewhat from the values vo found 
from the fit by Eq. ( fTSl i. This discrepancy may be due to the men- 
tioned observational errors, but it may also reflect an influence of 
effects neglected in the simple models given by the broadening 
functions A°.^ and A^.^. For instance the gradients of pulsational 
velocities or the stellar wind overimposed on the pulsations may 
contribute a distortion resembling P-Cyg shape to the line pro- 
file. More precise measurements as well as models of line for- 
mation will be needed to solve this question. 

In Table |2] we give the amplitude vi of the first periodic term 
in the expansion given by Eq. ( fT2b . It can be seen that this ampli- 
tude is smaller for the lines of ionized elements and it is slightly 
correlated with the equivalent width of the lines of neutral el- 
ements. This result indicates that the amplitude of pulsational 
velocities depends on the depth in the atmosphere where the 
line is formed (cf. Butler et al. 1996, Petterson et al. 2005, and 
references therein). This also qualitatively agrees with the cor- 
relation between line-depths and amplitudes of radial-velocity 
curves found in Cepheids by Nardetto et al. (2007 and 2009). 

Similarly, the line-strength factors (see Fig. [3] and on-line 
Table |4]i of the lines were expanded into the Fourier series 
and the amplitudes si of the first component are also given in 
Table |2] The negative values of this amplitude are assigned to 
the lines which are approximately in an antiphase with the other 
lines. These are the two lines of ionized elements in our set, for 
which (^1 ^ -35°, while for the lines of neutral atoms we have 
(pi ^ 165°. This amplitude (anti-)correlates well with the ex- 
citation potential e (+ ionization potential cf. Moore et al. 
1966) of the lower level of the line. The largest deviation can be 
seen for the Ca I line, which, however, may be blended with the 
Fe I (169) fine 6344.15 A with e = 2.42 eV. These variations of 
line strengths measured by the method of relative line photome- 
try (Hadrava 1997) thus enable one to easily find the changes of 
atmosphere temperature in the course of the Cepheid pulsation 
period (cf. Krockenberger et al., 1998, Kovtyukh and Gorlova 
2000). 

As a next step we performed disentangling of a wide spectral 
region (6506-6618 A) around the Ha line. We chose sampling 
in 1024 bins with a step of 5km/s per bin. The telluric lines were 
disentangled using a template obtained from spectra of 68u Her 
again. The pulsational velocities and line-strength factors ob- 
tained by disentangling of the whole spectral region into one 
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Fig. 4. Phase-dependence of pulsational velocity in km/s for 
spectral lines in the Ha region. The one-component solution is 
marked by (green) full circles, low- and high- excitation compo- 
nents of the two-component disentangling are marked by (blue) 
open circles and (red) crosses, resp. 
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Fig. 5. Phase-dependence of line-strength factors s for spectral 
lines in the Ha region. (The meaning of symbols is the same as 
in Fig. El) 

stellar component are given in Tables [3] and |4] (in columns la- 
beled H„3) and drawn (by full circles in green) in Fig. |4]and|5] 
respectively. The semiamplitude of the first Fourier mode of the 
pulsational velocity is vi - 21.40 km/s; the line-strength factors 
have a semiamplitude si = 0.062 and phase shift (0i ^ -20°) 
comparable to that of the ionized atoms in Fig. [3] The phase de- 
pendence of residual spectra is drawn in Fig.|6] (The wavelength 
A is drawn on a logarithmic scale and labeled in A.) It can be 
seen here that the strength of the Ha wide wings is almost in an- 
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6506 A 6618 



Fig. 6. Residual spectra in the Ha region after 1 -component (+ 
telluric spectrum) disentangling. The vertical offset is propor- 
tional to the phase. 



tiphase with its narrow core as well as with strengths of weaker 
metallic lines. The residuals in the wings are thus in absorption 
and the narrow lines in emission around phase 0.0 around max- 
imum expansion (cf. the bottom and uppermost lines) and it is 
reversed around phase 0.5 at the infall (the lines around the mid- 
dle of the Figure). Residuals of some of the narrow lines have a 
shape of P-Cyg or inverse P-Cyg profiles in some phases. This 
indicates that the Doppler shifts of these lines differ from the 
other lines. 

To decrease the residual spectrum, we must allow each line 
to vary its line strengths in agreement with the phase-locked 
changes of the atmosphere temperature. We thus also disentan- 
gled the Ha region into two independent components, starting 
with the two solutions from Table |2] which are the opposite ex- 
tremes in si, i.e. the solution for the line Fe I (13) 6358.69 A and 
Si II (2) 6347.10 A. The spectrum splits into two components, as 
it is shown by the upper two lines in Fig.|2] The third very bottom 
line in this Figure gives the telluric spectrum. It can be seen here 
that the Ho' line is contained mostly in the second component ob- 
tained from the initial approximation corresponding to the Si II 
line, while the narrow metallic lines are distributed between both 
components, each one to a different proportion. A small contri- 
bution to the core of Ha also appears in the first spectrum. This 
contribution has an inverse P-Cyg shape. 

We let the velocities and line-strengths of both these compo- 
nents in all spectra converge to the best fit. The (O-C)^ of resid- 
ual spectra decreased to 12% in comparison with the previous 
one-component (+ 1 telluric component) solution. The resulting 
velocities are labeled as H„i or H„2 in Tableland drawn in 
Fig. |4]by (blue) open circles or (red) crosses for the components 
developed from the solutions for neutral and ionized atoms, re- 
spectively. The first Fourier modes of the pulsational velocities 
are vi = 20.18 km/s or vi = 17.38 km/s, resp. The phase shift 
of the first mode of velocities converged to 0i ^ 228° for the 



II 

— I 1 1 1 ^ 1 1 1 1 1 1 

6550 ^ 6600 

Fig. 7. Spectrum of 6 Cep in the H,, region disentangled 
into two components. The top (blue) line corresponds to the 
low-excitation component, the middle (red) line to the high- 
excitation component and the bottom (black) line to the telluric 
component. 

Table 5. List of spectral lines disentangled in the Ha region 



^[A] 


multiplet 


X + e [eV] 


EW, [A] 


EW2 [A] 


6516.08 


Fe II (40) 


7.87 -1- 2.89 


0.046 


0.102 


6518.37 


Fe I (342) 


2.83 


0.069 


0.020 


6546.25 


Fe I (268) 


2.76 


0.087 


0.063 


6587.63 


C I (22) 


8.58 


-.004 


0.033 


6592.93 


Fe I (268) 


2.73 


0.097 


0.082 


6604.60 


Sell (19) 


6.54+1.36 


0.058 


0.064 



first component, unlike the values (pi ^ 197° + 3° found for the 
second component of Ho- region, its one-component solution, 
as well as for the solutions of individual lines listed in Table |2] 
(with the exception of the value ^ 203° for the line Fe II). The 
line-strength factors (given in Table |4] and Fig. |5]l have semi- 
amplitudes si = 0.577 with phase shift (0i - 166°) compa- 
rable to that of neutral atoms in Fig. [3] for the first component 
and si - -0.152 for the second component. The negative sign 
here denotes again the approximate antiphase {<pi ^ -20°, cf. 
Fig. Is). Let us note that the line C I (22) entirely appears and 
the lines Fe II (40) and Sc II (19) are more pronounced in the 
second component also containing the major part of the high ex- 
citation Ha line (e = 10.20 eV), while the other low excitation 
Fe I lines identified in Fig. [T] and Table |5] (where estimates of 
equivalent widths in the first and second component are given) 
are by greater or equal part contained in the first "low excitation" 
component of the spectrum. 

These results are consistent with those for individual lines 
presented above and with the variability of radial velocities and 
equivalent widths found using classical methods by different au- 
thors cited there. 



4. Conclusions 

Our measurements of pulsational velocities of S Cep using a 
newly generalized version KOREL09 of the code for Fourier 
disentangling and line-strength photometry confirmed that the 
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phase variations of velocities and strengths of different spectral 
lines depend on the depth of their formation in the stellar atmo- 
sphere and on the excitation potential of their lower levels. This 
means that any elaboration of the Baade - Wesselink method to a 
higher precision must take into account the structure and dynam- 
ics of Cepheid atmospheres, where the instantaneous effective 
stellar radii corresponding to either the photometric or interfer- 
ometric version of the method do not precisely follow the true 
motion of any particular layer of the atmosphere or velocities 
found from Doppler distortions of any particular stellar line. Our 
method of disentangling spectra provides a tool for observational 
probing of the structure of the pulsating stellar atmospheres. 

Originally designed and well proved for disentangling of 
spectra of multiple stars, our method can also be used for studies 
of Cepheids in binaries, which can yield more information about 
the basic physical parameters of the stars and thus also about 
the PL-relation in different conditions. Further sophistications 
of the method are possible, e.g. directly fitting the coefficients 
of the Fourier series given by Eq. ( fTSl i instead of independent 
velocities in individual exposures to distinguish the orbital and 
pulsational changes of radial velocities. These coefficients could 
be then used for the distance determination in the way introduced 
by Krockenberger et al. (1997) or for the classification according 
to Deb and Singh (2009). 

Another possible future improvement is to take into account 
for the pulsational broadening functions the results of radia- 
tive transfer and line formation in differentially pulsating atmo- 
spheres. This is particularly challenging for non-radial pulsa- 
tions, where the approximation of an unperturbed moving atmo- 
sphere is even less physically substantiated but even more com- 
monly used to avoid the complexity of the problem. Note that the 
output of residual spectra in rest wavelengths of any component 
or the centre of mass of a multiple system was implemented in 
the KOREL code by its author as a "first aid" for studying pul- 
sational or other perturbations of spectra not directly included 
in the model of the broadening functions A. However, such an 
approach should not be used as a black box without understand- 
ing its underlying assumptions, namely that the free parameters 
of the neglected effect (e.g. the pulsations) do not correlate with 
the parameters taken into account (e.g. orbital parameters). If 
we find a solution with a simplified model, it is not proof that 
its assumptions are correct, because neglecting a problem is not 
a true solution of the problem. The disentangling of pulsations 
in spectra of single or multiple stars thus requires one to include 
the pulsational variations in the model by which we fit the ob- 
servations. We have demonstrated such an approach here in its 
simplest form. A future sophistication will be needed to account 
for more accurate observations. 
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Table 3. Pulsational velocities v for individual lines (cf. Table |2| and Ha region at all exposures. 



HJD 

-2454000 



phase 6265.14 6343.71 6347.10 6355.04 6358.69 6400.01 6456.39 6663.45 H,,i H„2 H„ 



698.3027 .2426 -3.472 -4.805 -3.114 -4.231 -3.812 -3.377 -5.866 -4.813 -10.724 -3.538 -3.581 

702.3672 .0000 -19.120 -20.131 -18.591 -18.690 -18.623 -21.070 -17.497 -20.683 -20.926 -18.053 -21.937 

708.4863 .1403 -10.826 -10.841 -9.845 -10.457 -9.998 -10.685 -10.559 -11.224 -16.223 -9.298 -10.475 

708.4925 .1414 -10.311 -11.554 -9.525 -11.006 -10.281 -9.949 -9.568 -10.841 -16.093 -8.544 -10.079 

709.3067 .2932 .463 .161 .131 .168 .355 1.027 -1.491 .367 -6.903 -.535 .282 

709.3918 .3090 .541 1.024 .803 .661 1.017 2.076 -.660 1.090 -5.914 .429 1.306 

709.4934 .3279 3.044 3.655 2.640 2.909 2.310 3.953 1.411 3.139 -3.749 2.340 3.406 

710.2962 .4775 11.764 13.170 12.519 12.540 11.533 14.355 10.445 12.580 6.654 10.167 13.655 

710.3885 .4947 12.200 14.001 12.923 13.476 11.755 15.001 11.141 13.225 7.582 10.522 14.431 

713.3641 .0492 -16.299 -18.433 -15.865 -17.112 -14.856 -17.865 -15.884 -17.800 -21.821 -14.746 -17.823 

719.2701 .1498 -9.941 -10.207 -9.739 -9.979 -9.482 -10.167 -9.862 -10.687 -16.053 -9.057 -9.976 

719.3365 .1622 -8.472 -9.176 -8.075 -8.449 -7.727 -8.343 -8.468 -9.273 -14.641 -7.795 -8.541 

719.3823 .1707 -8.104 -6.857 -7.089 -6.962 -7.039 -7.396 -7.464 -8.359 -13.749 -7.148 -7.656 

719.4983 .1923 -6.140 -5.839 -5.920 -5.720 -5.739 -5.189 -5.965 -7.151 -12.214 -5.771 -5.955 

719.5025 .1931 -7.041 -6.427 -6.718 -6.484 -6.072 -6.175 -6.978 -7.242 -13.088 -6.391 -6.643 

719.5764 .2069 -5.763 -6.034 -4.739 -4.705 -4.399 -4.846 -6.002 -5.864 -11.677 -5.091 -5.187 

736.4506 .3513 5.086 4.986 4.263 4.868 4.232 5.656 2.834 4.452 -2.272 4.177 5.114 

737.2522 .5007 14.810 14.151 12.956 13.590 12.394 15.191 11.347 15.769 9.054 11.995 15.898 

737.4991 .5467 15.897 16.675 15.857 15.609 14.873 17.766 13.501 16.647 11.422 13.820 18.197 

737.6202 .5693 18.153 18.430 17.510 17.579 16.354 19.509 15.891 18.914 13.497 15.659 20.182 

738.2559 .6877 25.667 26.431 25.142 26.082 24.013 26.733 24.830 26.803 21.757 23.688 27.721 

738.3353 .7025 25.542 26.555 24.756 25.440 24.145 26.485 25.407 27.217 23.200 24.194 28.507 

738.4407 .7222 22.621 24.298 22.876 23.889 22.326 24.393 24.408 23.966 23.626 22.297 27.316 

738.5058 .7343 21.891 22.284 22.006 22.177 21.452 22.448 23.506 23.296 24.665 20.887 26.895 

738.5121 .7355 22.000 22.983 21.673 22.090 21.350 22.549 23.257 23.523 24.795 20.794 26.921 

738.5852 .7491 18.959 19.164 18.943 19.028 18.171 18.670 20.746 20.229 24.993 18.405 24.877 

744.3990 .8325 -12.053 -12.730 -12.948 -13.166 -10.740 -14.169 -10.648 -14.182 15.896 -11.724 -16.115 

751.2407 .1074 -13.394 -13.437 -12.515 -12.752 -12.051 -13.271 -11.996 -13.347 -17.913 -11.563 -13.046 

751.2445 .1081 -13.639 -12.864 -11.994 -12.394 -11.397 -12.946 -11.251 -12.639 -17.725 -11.017 -12.618 

752.2419 .2940 -.351 -.150 -.213 -.023 -.118 .540 -1.416 .192 -7.474 .391 .310 

752.3239 .3093 .104 .517 .161 .235 .368 1.205 -1.246 .242 -6.694 .390 .755 

752.3612 .3162 .601 .873 .901 .883 1.079 1.952 -.628 .618 -5.854 .430 1.406 

757.2615 .2294 -4.640 -4.563 -4.879 -4.458 -4.136 -4.008 -5.363 -4.740 -11.266 -4.572 -4.490 

758.3017 .4232 8.058 8.891 8.077 8.648 7.927 9.663 6.442 8.495 2.435 6.836 9.369 

761.4677 .0132 -19.167 -19.361 -17.262 -17.943 -17.526 -19.989 -15.940 -19.126 -22.674 -15.931 -20.340 

771.2950 .8445 -16.757 -18.029 -17.489 -17.553 -14.807 -20.774 -15.062 -18.224 10.474 -15.843 -21.765 

771.3017 .8457 -17.270 -17.912 -17.540 -18.093 -17.744 -20.982 -14.826 -18.092 9.583 -16.326 -22.067 

772.2613 .0245 -18.332 -18.172 -17.079 -18.665 -17.204 -19.081 -15.587 -18.398 -21.602 -15.372 -19.205 

774.2135 .3883 6.001 5.869 5.153 5.797 5.404 7.304 4.523 6.392 -.147 5.395 7.202 

776.1982 .7582 14.591 15.493 14.505 14.258 14.348 13.684 16.629 15.125 23.436 14.231 21.215 

776.2430 .7665 11.360 11.644 11.325 11.487 10.499 10.375 13.564 12.644 23.432 11.784 18.674 

776.2868 .7747 7.889 7.464 7.835 7.097 7.709 6.312 10.046 8.428 22.699 8.657 14.885 

776.2972 .7766 7.360 6.757 7.886 7.655 7.222 6.257 9.713 8.030 23.566 8.493 14.584 

781.2311 .6960 24.209 24.341 21.721 23.887 22.540 24.820 23.302 25.921 21.316 22.742 26.844 

798.3846 .8925 -22.209 -20.912 -22.504 -22.707 -24.052 -26.179 -20.384 -24.487 -4.586 -21.204 -28.487 

798.3939 .8943 -22.432 -21.594 -21.541 -21.999 -20.465 -25.455 -19.844 -23.852 -4.747 -20.621 -27.909 

798.4495 .9046 -22.934 -23.639 -22.419 -23.410 -20.900 -26.265 -20.292 -24.475 -7.728 -21.228 -28.217 

824.3317 .7277 18.792 20.530 17.753 18.852 17.995 19.383 19.651 20.842 20.695 16.979 23.366 

843.2488 .2528 -2.825 -2.918 -3.969 -3.555 -2.882 -2.933 -4.288 -2.926 -10.146 -2.626 -2.791 

857.2159 .8555 -20.871 -20.889 -20.036 -20.335 -20.987 -23.105 -17.185 -20.980 5.365 -19.282 -25.892 

910.6409 .8111 -7.020 -7.160 -6.173 -6.396 -5.840 -9.023 -5.008 -7.355 20.318 -5.737 -6.244 

923.5714 .2207 -5.783 -5.871 -5.055 -5.467 -5.484 -5.258 -6.204 -5.840 -12.770 -5.549 -5.808 

923.6211 .2299 -5.610 -5.789 -4.929 -5.372 -5.455 -4.547 -6.149 -5.344 -12.815 -5.218 -5.527 

927.5339 .9591 -23.305 -22.027 -21.087 -23.797 -19.956 -24.735 -19.891 -23.971 -21.156 -19.528 -25.434 

928.4066 .1217 -12.945 -13.761 -11.411 -12.335 -12.105 -13.326 -12.596 -12.991 -18.469 -10.608 -12.469 

928.4523 .1302 -12.809 -13.317 -10.827 -11.937 -11.876 -12.539 -11.972 -13.189 -18.086 -10.433 -12.017 

928.4885 .1370 -12.083 -12.143 -10.498 -11.798 -11.156 -11.789 -11.150 -12.400 -17.737 -10.201 -11.629 

928.6318 .1637 -9.680 -9.726 -9.079 -9.399 -9.103 -9.293 -9.046 -9.384 -14.967 -8.243 -8.868 

931.4636 .6914 25.727 25.717 24.632 26.043 24.502 26.672 25.196 26.687 21.251 23.321 27.381 

931.6042 .7176 24.333 24.484 22.787 24.512 22.557 24.579 24.155 25.221 23.197 22.469 27.464 

931.6139 .7194 23.349 23.641 22.779 23.520 22.059 23.678 23.775 24.484 23.100 21.896 27.099 

932.5844 .9002 -22.441 -22.897 -22.678 -23.743 -20.409 -25.682 -19.925 -24.120 -7.148 -20.443 -27.991 

932.6097 .9049 -23.106 -24.767 -22.705 -23.376 -20.218 -25.918 -19.832 -24.021 -8.526 -20.536 -27.756 

932.6138 .9057 -23.617 -25.700 -22.690 -23.541 -22.236 -26.092 -19.773 -23.632 -8.680 -20.248 -27.686 

932.6178 .9064 -21.966 -24.284 -22.814 -23.471 -19.347 -26.057 -20.040 -24.018 -8.634 -20.353 -27.959 

933.5557 .0812 -15.562 -15.682 -15.102 -15.644 -14.291 -16.345 -14.030 -15.348 -20.031 -12.778 -15.282 

934.5024 .2576 -3.043 -3.294 -2.396 -2.486 -3.021 -2.173 -4.223 -2.895 -10.500 -2.929 -2.991 

937.5199 .8199 -10.328 -12.313 -10.500 -10.560 -10.105 -13.348 -8.644 -11.404 16.376 -9.508 -11.648 
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Table 4. Line-strength factors for individual lines (cf. Table|2]i and Ha region at all exposures. 
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